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Abstract 

A new theoretical framework is presented that uses a localized triaxially deformed Gaussian as the 
single-particle wave packet. The model space enables us to adequately describe both the deformed 
mean-field structure and the cluster structure within the same framework. The improvement over 
the original version of antisymmetrized molecular dynamics which uses the spherical Gaussian is 
verified by its application to the 20 Ne nucleus. The coexistence and interplay of the cluster structure 
and the deformed mean-field structure in the low-excited states of 20 Ne is studied. In particular, 
the intra-band E2 transition probabilities in K w =0f and 2~ bands are reproduced without any 
effective charge. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



In light nuclei up to the beginning of the s<i-shell, the cluster structure is one of the most 
essential features of the nuclear structure together with a mean-field-like (shell-model-like) 
structure. On the heels of the cluster structure studies of the p-sd shell nuclei many 
theoretical studies have been conducted to investigate the cluster aspects of the heavier 
«*,/ shell nuclei QQ. Experimentally, some negativity levels which have large a- 
spectroscopic factors have been observed in 40 Ca and 44 Ti jj], and these states have been 
considered to be K n =0~ band members of the parity doublet bands with K n =0 + which have 
an a+core structure. It is of great interest and importance to survey the cluster aspect of 
these heavier isotopes, but the study seems to be in a rather intermediate stage. Another new 
field for cluster study is the unstable nucleus. Due to advances in experimental technique, 
many new aspects and phenomena such as neutron halo, neutron skin and the change of 
the magic number P, |^] have been observed and studied, and now an increasing body of 
information on the unstable nuclei is available. It has been shown by many theoretical and 
experimental studies that the cluster-like correlation is also important and the exotic cluster 
structure which cannot be seen in stable isotopes appears in light neutron-rich isotopes. 
Therefore, new cluster aspects are expected in heavier stable and unstable nuclei. 

One of the difficulties of cluster structure study in heavier nuclei is the existence and 
importance of the mean-field deformation of the nucleus, which should be properly treated 
together with clustering. It is quite conceivable that the mean-field deformation and cluster- 
ing coexist, compete or are mixed. For the study of the nucleus in such situations, we have 
to treat the cluster structure and the deformed mean-field structure within the same frame- 
work. Another difficulty is the increasing importance of the spin-orbit force. The strong 
effect of the spin-orbit force will dissolve the cluster structure. Therefore the expectation 
value of the spin-orbit force can be regarded as a kind of measure that tells whether or 
not a cluster structure with doubly closed shell subunits develops. The importance of the 
spin-orbit force increases as we proceed to heavier nuclei. 

To study the cluster aspects of heavier nuclei, we have introduced deformed-basis antisym- 
metrized molecular dynamics (deformed-basis AMD), which employs a localized triaxially 
deformed Gaussian as the single-particle wave packet j^. In this work, we present the de- 
tailed framework of deformed-basis AMD and argue the advantages of this new framework. 
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In short, the present framework has two advantages. The first is its better description of the 
deformed mean-field achieved by making the single-particle wave packets deformable. The 
second is its improved incorporation of spin-orbit force effects. These advantages become 
clear in its application to the low-lying rotational bands of the 20 Ne nucleus. 

The K n =0f bands of 20 Ne are regarded as parity doublet bands with an a+ O cluster 
structure. However, at the same time, the impurity of the a+ 16 cluster structure in the 
ground band has been discussed by many authors, while the K n =0~ band is considered 
to have an almost pure a+ 16 cluster structure. Indeed the intra-band B(E2) values of 
the ground band have been underestimated by a+ 16 cluster models. It is considered 
that this deficiency originates in the mixture of the deformed mean-field structure with an 
a+ 16 cluster structure. The lowest negative-parity band is the K 7T =2~ band which has a 
(0p)~ 1 (sd) 5 structure. Unlike K 7r =Q^ 1 bands, this band should be understood in terms of 
the deformed mean-field. Thus, in the low-lying states of 20 Ne, the cluster structure and 
the deformed mean-field structure coexist and are mixed with each other. This intermingled 
situation has not been treated in a consistent manner within a full microscopic framework. 
It will be shown that the present framework is flexible enough to describe the coexistence 
and the mixture of the cluster structure and the deformed mean-field structure. 

The contents of this article are as follows. In the next section, the framework of deformed- 
basis AMD is given and the calculational procedure adopted in this study is explained. The 
features and the advantages of this new framework are explained in the section IIHI The 
advantages in the practical calculation are shown by applying the method to the 20 Ne nucleus 
(section HVjl . In the course of this application, the characters of the low-lying rotational 
bands of 20 Ne, K n =0f,0^,0^ and 2" bands are discussed, and the calculated results are 
examined. In the last section, we summarize this work. 



II. FRAMEWORK OF DEFORMED-BASIS AMD 

In this section, the framework of deformed-basis AMD is presented and the calculational 
procedure adopted in this study is explained. In the following, we call usual AMD, which 
uses a spherical Gaussian as the single-particle wave packet 'spherical-basis AMD' to keep 



the distinction clear. For a more detai 
the reader is directed to references 



explanation of spherical-basis AMD framework, 

mm 
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A. Wave function and Hamiltonian 



In deformed-basis AMD, the intrinsic wave function of the system with mass A is given 
by a Slater determinant of single-particle wave packets; 

®int = -j==A{(pi,ip 2 ,...,(PA}, (1) 
<Pi(r) = </>i(r)xi£i, (2) 

where <fi is the ith single-particle wave packet consisting of spatial 0«, spin x% an d isospin 
£i parts. Deformed-basis AMD employs the triaxially deformed Gaussian centered at Zj 
as the spatial part of the single-particle wave packet, while spherical-basis AMD limits the 
Gaussian to spherical shape: 



eXP ] - Z ^) 2 p 



a=x,y,z 

|2 i 1/3 12 



Xi = atiXi + PiXl* H + \Pi\ 
£i = proton or neutron. (3) 

Here, the complex number parameter Zj which represents the center of the Gaussian in phase 
space takes an independent value for each nucleon. The width parameters u x , v y and v z are 
real number parameters and take independent values for each direction, but are common to 
all nucleons. Spin part Xi is parametrized by and $ and isospin part & is fixed to up 
(proton) or down (neutron). Zj, v x ,v y , v z and ccj, /?» are the variational parameters and are 
optimized by the method of frictional cooling explained in the next subsection. By using 
a deformed Gaussian basis, we can improve the description of deformed nuclei. The effects 
resulting from the deformed Gaussian basis are explained in the section IIHI and are shown 
in its practical application to the 20 Ne (section llVjl . As the variational wave function, we 
employ the parity projected wave function in the same way as many other AMD studies 

$ ± = P^int = i ^^^nU (4) 



here P x is the parity operator and is the intrinsic wave function given in EqQ. 
The Hamiltonian used in this study is as follows; 

H = f + V n + V c -f g , (5) 



where T and T g are the kinetic energy and the energy of the center of mass motion, respec- 
tively. The expectation value of the T g is evaluated without any approximation, but the 
wave function contains the spurious component of the center-of-mass motion. Because the 
wave function of the center-of-mass also has a deformed Gaussian form, it becomes insep- 
arable when the angular momentum is projected and/or the multiple Slater determinants 
are superposed, prescription should be developed. In this study, we have approximated 
its influence by subtracting the expectation value of the center-of-mass kinetic energy. We 
have used the Gogny force D1S force as an effective nuclear force V n . Coulomb force V c is 
approximated by the sum of seven Gaussians. 



B. Energy variation, angular momentum projection and generator coordinate 
method 

We perform the variational calculation and optimize the variational parameters included 
in the trial wave function Eq(j3J) to find the state that minimizes the energy of the system 

^ = ^L/, H = H + Vm * (6) 

We add the constraint potential V cns tto the Hamiltonian H to obtain the minimum energy 
state under the optional constraint condition. In this study, we constrain matter quadrupole 
deformation by employing the potential V cnst = v cnst ((/3 2 — (3q) 2 and we obtain the optimized 
wave function $ ± ( / 9o) = P^^intiPo) as a function of the deformation parameter (3 . The 
evaluation of the quadrupole deformation parameter (3 is explained in reference . At the 
end of the variational calculation, the expectation value of V cns t should be zero in principle, 
and in the practical calculations we confirm that it is less than 100 eV. 

Energy variation of the AMD wave function is performed by the frictional cooling method, 
which is one of the imaginary time development methods. The reader is directed to references 
[l^l Q| for a more detailed description. The time development equation for the complex 
number parameters Zj, and $ is as follows; 

dX t _^ d (*±\H\*±) 

dt ~ hdx* <$±|$±> ' ^--^'-'^ W 
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and for the real number parameters v x , v y and v z \ 

dv a fjf d ($ ± |H|$ ± ) 



(a = x,y,z) (8) 



dt hdu a <$±|$±) 

Here X$ is Zj, oii, or fy. \x and // are arbitrary negative real numbers. It is easy to show 
that the energy of the system decreases as time develops, and after sufficient time steps we 
obtain the minimum energy state. 

From the optimized wave function, we project out the eigenstate of the total angular 
momentum J, 

*mM) = Pmk&W = Pm^UPo). (9) 

Here Pmk ^ s the total angular momentum projector. The integrals over the three Euler 
angles included in the P(, IK are evaluated by numerical integration. 

Furthermore, we superpose the wave functions & J A f K which have the same parity and 
angular momentum but have different values of deformation parameters (3q and K. Thus 
the final wave function of the system becomes as follows: 

where quantum numbers other than total angular momentum and parity are represented by 
n. The coefficients c n , c' n ,... are determined by the Hill- Wheeler equation, 

5(« ± I^K ± )-en« ± K ± ))=0. (11) 

It is sometime useful to include the intrinsic wave function obtained from the variational 
calculation after the projection to the opposite parity. For instance, suppose that $i nt (+)((3o) 
and $j„t(_)(/3o) are the intrinsic wave functions obtained by the variational calculation after 
the projection to positive- and negative-parity, respectively. We include P^K^inti-) (A)) i n 
the GCM basis as well as PMK^int(+)(Po) to describe the positive-parity state $^ + ; 

$f = c n Pj£ K <$> int(+) ((3 Q ) + c' n P^ K ^ int{+) ((3' ) + ■■■ 

+4P£ K »* M (- ) (PS) + CnPMK^nt^) + ' • • ■ (12) 

In the present study, we have employed all of the intrinsic wave functions on the obtained 
positive- and negative-parity energy curves as the GCM basis. 
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III. FEATURES AND ADVANTAGES OF THE PRESENT MODEL 



The framework of AMD has some advantages over other theoretical methods. For in- 
stance, parity projection before the variation, angular momentum projection after (before) 
the variation and superposition of the Slater determinants (GCM) are easily executable 
in this framework. It is well known that parity projection is essentially important to de- 
scribe the parity asymmetric shape of the nucleus and parity asymmetric cluster structure 
14, 3|. Angular momentum projection and superposition of Slater determinants are indis- 



pensable to study the level scheme of the nucleus and the nuclear-structure-change in each 
state beyond the mean-field picture. In this section, besides these advantages which are 
common to spherical- and deformed-basis AMD, we explain the features and advantages of 
deformed-basis AMD in the study of the heavier nuclei. First, we discuss the description of 
the deformed mean-field structure and cluster structure in this model. Then we discuss the 
evaluation of the spin-orbit force. 



A. Deformed mean-field structure and cluster structure 

In deformed-basis AMD, deformed mean-field structure is described by deformed single- 
particle wave packets and the centroids of single-particle wave packet Zj gathered to the 
center of the nucleus, while the cluster structure is described by the localization of Z^s into 
several parts of the coordinate space. 

Roughly speaking, in deformed-basis AMD, the change of nuclear shape as a function of 
nuclear deformation is classified into three patterns as sketched in FIG.^ In this figure, for 
the sake of simplicity, we consider only the elliptical deformation of the mean-field and the 
cluster structure consisting of two spherical subunits. The first pattern shows deformation 
of the mean- field (FIG. (i)). When the nuclear deformation becomes larger, only the 
deformation of the mean-field becomes larger. In this case, all Zj in EqfjSJ) are gathered 
in the small region around the center of nucleus and the single-particle wave packets are 
deformed. In this limit, the nuclear deformation is described by the deformation of the 
single-particle wave packets and the strong effect of the Pauli principle (single particle motion 
in the deformed mean- field). In other words, the intrinsic wave function Eq. ([T} contains 
a large amount of the deformed shell model configurations. This structure is not described 
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>► large deformation 



FIG. 1: Schematic figure showing three different patterns of the nuclear deformation in deformed- 
basis AMD. Points in the figure indicate the centroids of the single-particle wave packets (Re(Zj)) 
of a deformed-basis AMD wave function. 

well by spherical-basis AMD and this is the reason that spherical-basis AMD sometimes 
provides an inadequate description of deformed heavier nuclei ^(|. The second pattern 
(FIG.[n (ii)) shows the deformation caused by clustering. In this pattern, the nucleus splits 
into two clusters. This structure is described by the localization of Zj into two parts and 
the spherical shape of the single-particle wave packets. As the nuclear deformation becomes 
larger, the distance between two clusters becomes larger, but the single-particle wave packets 
remain spherical. The third pattern (FIG. Q (hi)) has a mixed character of (i) and (ii). At 
small deformations, the nucleus has a mean-field character, but as the deformation becomes 
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larger, a cluster-like structure grows and the nucleus has an intermediate character. At large 
deformations, it changes to an almost pure cluster structure. In the application to 20 Ne, 
it will be shown that J 7r =2~(A' 7r =2~),l~(A' 7r =0~) and + (A' 7r =0 + ) curves correspond to 
patterns (i),(ii) and (iii), respectively. 

B. Evaluation of the spin-orbit force 

In medium-heavy nuclei, the importance of the spin-orbit force increases. In short, the 
spin-orbit force acts to dissolve the cluster structure. Since spherical-basis AMD does not 
assume any cluster structure, it can describe the dissociation of the cluster structure due 
to the spin-orbit force. However, it was pointed out that in heavier isotopes, the spherical- 
basis AMD wave function provides a much smaller expectation value of the spin-orbit force 
compared to HF theory when the same interaction Skyrme III force is adopted Q|. We 
consider that this is because single-particle wave packets are limited to the spherical Gaussian 
form in spherical-basis AMD. Since the spin-orbit force contains the derivative operator, it 
is sensitive to the slope and the deformation of the single-particle wave packet. When the 
system favors the mean- field deformation, doubly closed-shell clusters such as a and 16 O 
clusters are difficult to form. This often means that deformed-basis AMD gives us larger 
expectation values of the spin-orbit force than spherical-basis AMD does, because the latter 
tends to give us a larger amount of the clustering component in the system wave function 
than the former. For instance, deformed-basis AMD gives twice as large an expectation value 
of the spin-orbit force (-19.0 MeV) as spherical-basis AMD gives (-11.6 MeV) in the 24 Mg 
nucleus. In contrast, when a nucleus favors a cluster structure consisting of doubly-closed- 
shell subunits like A 7r =0~ band of 20 Ne, spin-orbit force does not act at all and the shape 
of the single-particle wave packet is not distorted by it. Therefore, the expectation value of 
the spin-orbit force can be regarded as a type of measure that identifies the formation and 
dissociation of the cluster structure. In the next section, it will be shown that the different 
characters of the rotational bands of 20 Ne are directly reflected in the expectation value of 
the spin-orbit force. 
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IV. APPLICATION TO 20 Ne AND DISCUSSION 



We have applied the present framework to the 20 Ne nucleus to study the coexistence and 
mixing of the cluster structure and the deformed mean-field structure. First, we investigate 
the structure change of the wave functions as a function of matter quadrupole deformation. 
Then we study the detailed properties of each band. 

A. Overview of the energy curves 

After the variational calculation and the angular momentum projection to the J 7T =0 + 
(K n =0 + ), 1~ (K n =0~) and 2~ (K 7T =2~), we obtain three energy curves as functions of the 
deformation parameter (3 (FIG. |2J). In the positive-parity state, we obtain the J n =0 + curve 
(solid line with circles in FIG. |2J), and when we perform the GCM calculation, J 7r =0 1 l " and 
J T =04 states are mainly obtained from this curve. Similarly, J' K =2\ and J 1T =2~l states are 
obtained by superposing the wave functions on the J 7r =2 + (K n =0 + ) curve and so on for 
other higher spin states. In the negative-parity state, we obtain the J w =l~ curve (dashed 
line with triangle) and 2~ curve (dotted line with boxes) which mainly contribute to the 
J*=l- (K n =0^) and J n =2- (K v =2~) states, respectively. Detailed results of the GCM 
calculation are given in the next subsection, and in this subsection, we concentrate on the 
J 7r =0 + , 0~ and 2~ curves to discuss the different nature of the wave functions. The results 
obtained by spherical-basis AMD are also plotted for comparison (lines without marks). 
Three patterns of the nuclear shape change mentioned in the previous section can be seen 
on these three energy curves obtained by deformed-basis AMD. When we compare the 
J w =l~ (K w =0~) curve of deformed-basis AMD with that of spherical-basis AMD, they are 
almost identical, while there are differences in the position of the minimum energy point 
and in the minimum energy in the J n =2~ (K 7T =2~) curves. In deformed-basis AMD, the 
J n =2~ curve has a minimum energy point around /?=0.45, but in spherical-basis AMD, it 
has an unclear minimum energy point around (3=0.25 and its energy is underestimated by 
about 5 MeV compared with deformed-basis AMD. These similarity and difference are due 
to the different nature of the J n =l~ and 2~ curves. The J 7r =l~ curve has an almost pure 
cluster structure, while the J n =2~~ curve has a deformed mean-field structure. Since the 
cluster structure is well described by both spherical- and deformed-basis AMD, there is no 
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FIG. 2: /3-energy curves of 20 Ne (energy curves as a function of the deformation parameter 0). 
Lines with symbols (Lines) are the result of deformed-basis AMD (spherical-basis AMD). K 7T =0 + 
and K 7T =2~ states are more deeply bound in deformed-basis AMD than in spherical-basis AMD 
(See text). 



difference in their description of the J n =l~ curve. However, because deformed-basis AMD 
can describe the deformed mean-field structure better than spherical-basis AMD, there is a 
difference in the description of the J 7r =2~ curve. The different nature of these two states 
becomes clear in the deformation of the single-particle wave packets and in the distributions 
of the centroids of the single-particle wave packets. In the J 7r =l~ curve (dashed line in FIG. 
0J), single-particle wave packets are always almost spherical but the centroids of the twenty 
single-particle wave packets (black points in FIG. 0] (d) and (e)) are widely separated into 
two sub units, a (four points degenerate on the upper side) and 16 O (sixteen points are almost 
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FIG. 3: The ratio of the width parameter of the single-particle wave packet (y Vzjvx) as a function 
of the deformation parameter /3. In the present calculation, the single-particle wave packets are 
always almost axially symmetric and z axis is taken to be the symmetry axis. 

degenerate in the middle). As the nuclear deformation becomes larger, only the distance 
between two subunits becomes larger. Therefore, this nuclear shape change corresponds to 
the second pattern in FIG. In contrast, in the J 7r =2~ curve, the single-particle wave 
packets are significantly deformed, as a function of the nuclear deformation (dotted line in 
FIG.|2J). All of the centroids of the twenty single-particle wave packets are gathered around 
the center of the nucleus (FIG.|H(f) and (g)), even though the nuclear deformation becomes 
larger. We note that though the nuclear deformation is larger in wave function (g) than in 
wave function (d), the centers of the single-particle wave packets are crowded into the center 
of the nucleus in wave function (g). Therefore, this state corresponds to the first pattern in 
FIG. Hand shows a deformed mean-field nature. These different characters of the 7^=1" 
and 2~ curves are directly reflected in the expectation value of the spin-orbit force (FIG. 03). 
The spin-orbit force does not act at all in this state (dashed line with triangle in FIG. EJ). 
However, it acts strongly in the 7^=2" wave function (dotted line with boxes) which has an 
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FIG. 4: The density distributions of the intrinsic wave functions on the energy curve obtained by 
the deformed-basis AMD. Twenty centroids of single-particle wave packets (Re Zj, i=l~20) are 
denoted by the black points in each panel, (a), (b) and (c) are for the J 7r =0 + (K 7T =0 + ) state, (d) 
and (e) are for the J n =l~ (K w =0~) states, and (f) and (g) are for the J=2~ (K n =2~) states. 



(0p)~ 1 (sc/) 5 structure. The (0p)~ 1 (s<i) 5 structure of the calculated J 7r =2~ state is confirmed 
by investigating the structure of the single-particle orbit. The technique for studying a 
single-particle orbit structure within the AMD framework was developed by Dote, et al (l7|. 

The energy curve for the J 7r =0 + state is more deeply bound about 5 MeV at the mini- 
mum energy point in deformed-basis AMD than in spherical-basis AMD. The intermediate 
character of the J 7r =0 + curve is confirmed in the deformation and distribution of single- 
particle wave packets. In this curve, the deformation of single-particle wave packets (solid 
line in FIG. 0} becomes larger as the nuclear deformation becomes larger in the small- and 
medium- deformation region {(3 < 0.5). In this stage, the centroids of the single-particle 
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FIG. 5: The expectation value of the spin-orbit force for each parity projected wave function 
$ are plotted as functions of the deformation parameter (5. Lines with symbols (Lines) are 

the result of the deformed-basis AMD (spherical-basis AMD). 

wave packets are gathered around the center of the nucleus (FIG. Ufa)) and the nucleus has 
a deformed mean-field character. But the deformation of the single-particle wave packets is 
saturated around (3 ~ 0.5. At this stage, the centroids of the single-particle wave packets 
are separated into a and 16 parts (FIG.Efb)), though the distance between them is rather 
small (about lfm). Therefore, the nucleus has a good amount of a+ 16 cluster component, 
although the cluster structure is distorted by the deformed mean-field structure. After this 
stage, the distance between the two subunits becomes larger, as the nuclear deformation 
becomes larger, while the deformation of the single-particle wave packets does not change. 
Then around (3 ~ 0.65, the distortion eases and the single-particle wave packets rapidly 
become spherical. Thus, the J w =0 + curve corresponds to the third pattern in FIG.^ This 
change of nuclear shape and structure as a function of the deformation parameter (3 takes 
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FIG. 6: The excitation energies of the low-lying states of 20 Ne. The observed values (the results 
of the deformed-basis AMD+GCM calculation) are given by the solid lines (dashed lines). The 
spin-parity of each state is also shown. 



place gradually in the case of deformed-basis AMD. Therefore, the expectation value of the 
spin-orbit force decreases uniformly in deformed-basis AMD, while in spherical-basis AMD 
it becomes zero rather rapidly. 



B. Low-lying level structure in Ne 

After the angular momentum projection, we perform the GCM calculation by superposing 
all of the wave functions on the obtained energy curves, and calculate the level scheme of 
20 Ne (FIG. We have obtained two K n =0 + bands in positive-parity and K n =2~ and 
K n =0± bands in negative-parity. In this subsection, we discuss the detailed character of 
these bands and show the related quantities in due order. 
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TABLE I: Observed and calculated a-RW multiplied by 100 at the channel radius a = 6fm, 
or K w =0f, O4" and 07/ band members. For comparison, the results of the (sd) shell mod el ( SM) 
3, a+ 16 RGM (RGM) Q and (a+ 16 0)+( 8 Be + 12 C) coupled channel OCM (OCM) 

shown. 
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TABLE II: The squared amplitude of the a+ 16 component W (see appendix ^ and the expec- 
tation value of the spin-orbit force (V/ s ) (in MeV) of the GCM wave function for each state. 
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6 4 + 


0.67 (0.37) 


-6.8 




sr 




-14.4 




8 4 + 


0.55 (0.38) 


-7.4 




er 




-16.5 



1(3 



1. K n = Of band 



Experimentally, the band is built upon the 1~ state at 5.78 MeV, that is 1.05 

MeV above the a+ 16 threshold. In the present calculation, the excitation energy of the 1~ 
state is 6.72 MeV and the excitation energies of the higher spin states show good agreement 
with the observed values. 

The K n =0i band is regarded as having an almost pure a+ 16 cluster structure, since 
the band members have large a decay widths comparable with the Wigner limit. Indeed, in 
the present calculation, the obtained a decay widths are large enough to be comparable with 
the observed values (Table. HJ) and the wave functions of the member states have quite large 
amounts of the a+ 16 component (Table. EJ). The definitions and the method to evaluate 
these quantities are given in the Appendix. The almost pure cluster structure of this band 
is due to a dominant contribution from the wave functions on the K w =0~ (J 7r =l~,3~, ...) 
energy curves. As is discussed in the previous subsection, the wave functions on the J n —0~ 
energy curve have an almost pure a+ 16 cluster structure (and this character is common 
to other J n curves with K n =0~). As a consequence, the member states of this band have 
an almost pure a+ 16 cluster structure and small expectation values of the spin-orbit force 
(Table. Though we did not make any assumptions about the cluster structure in our 
calculation, the band members have almost pure cluster structures. This fact means that 
the cluster structure can exist as an independent degree of freedom of nuclear excitation, 
even though the effects of the spin-orbit force and the deformation of the mean-field are 
taken into account in the model space. 

The relative motions between a and 16 O clusters were also investigated following the 
traditional cluster model studies. The a reduced width amplitudes 3^l(o) (a-RWA) are 
shown in Fig. [7| (dotted lines). The definitions and the calculational procedure for the a- 
RWA are also presented in the appendix B As in other cluster studies, it is clear that K n —0^ 
band members belong to the N = 2n + L = 9 band members of the a+ 16 cluster structure. 
Here, N denotes the principal quantum number of relative motion and n denotes the radial 
quantum number (number of nodes). We also see that relative motion is suppressed inside of 
the nucleus and enhanced outside, which induces large intra-band E2 transition probabilities 
(Table. UTIJ). 
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FIG. 7: The a-RWA for the member states of the K n =0^ (solid lines), (dashed lines) and X 
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(dotted lines) bands. See text and the appendix B 



TABLE III: Observed (EXP) and calculated (AMD) intra-band E2 transition probabilities 



B(£2;Jf -> JJ) [e 2 fm 4 
multicluster GCM (5a) 



within the K w =0^ band. For comparison, a- 16 RGM (RGM) jl<| and 
221 ] are also shown. The effective charges which are used in each calculation 



are given in the bottom line. 



TTT . TTT 

J i ~> J f 


EXP 


RGM 


5aGCM 


AMD 


3~ 1- 


164 ± 26 


121.4 


155 


151.2 


5- -> 3- 




133.399 


206 


182.4 


T -► 5~ 




122.74 




141.6 


9" -► 7" 




67.83 




87.9 


5e/e 














2. Ground band (K 7 " = 0+ band) 

The obtained binding energy of the ground state is -159.8 MeV, which underestimates the 
observed value by about 0.8 MeV. The rotational character of the 0+, 2f and 4f states and 
the deviation from the rotational spectra of the 6f and 8^ states are reproduced, though 
the obtained excitation energies are smaller than the observed values. 

The ground band must be the parity doublet partner of the K n =0i band, since the 
ground band is the only K 7T =0 + band below the K v ={)^ band. However, at the same time, 
the existence of the deformed mean-field component in the ground band has been discussed. 
For example, if the pure a+ 16 cluster structure of this band is assumed, the intra-band E2 
transition probabilities are underestimated. Shell model calculation has shown the existence 
of broken symmetry components in the wave function that activates the spin-orbit force. 

The GCM wave functions describing the member states of this band have a mixed char- 
acter embodying both the cluster structure and the deformed mean-field structure. They 
have a dominant contribution from the wave functions on the K 7T =0 + (J 7T =0 + , 2 + , ...) curves. 
The GCM-amplitude of the Of state is shown in FIG. |H] (solid line). The GCM-amplitude 
that serves as a kind of measurement of the contribution to the GCM wave function from 
each intrinsic wave function is denned as 

Ic^A,)! 2 = |<<ei*M*(A)))| 2 , (is) 

where S?^* is the GCM wave function (EqJTUJ) and &mk(Po) ^ s ^ ne P ar ity and angular mo- 
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TABLE IV: Observed (EXP) and calculated (AMD) intra-band E2 transition probabilities 
B(£2;Jf -> JJ) [e 2 fm 4 ] within the K"=gt band. For comparison, the results of the (sd) 4 shell 
model (SM) 



li, a+ 16 RGM (RGM) 



ISO, (a+ lb O)+( 8 Be+ 12 C) coupled channel OCM (OCM) 



2(3 ] a 

each calculation are given in the bottom line. 



TTT TTT 
J i -f J f 


EXP 


SM 


RGM 


OCM 


5aGCM 


AMD 


2+^0+ 


65 ±3 


57.0 


36.2 


57.0 


50.0 


70.3 


4+^2+ 


71 ±6 


69.9 


45.22 


70.9 


64.4 


83.7 


6+^4+ 


64 ± 10 


57.9 


36.5 


57.1 


55.3 


52.7 


8+^6+ 


29 ±4 


35.5 


19.7 


34.8 




21.0 


Se/e 




0.54 





0.155 









mentum projected wave function (EqJUj). The GCM-amplitude reaches a maximum around 
(3 ~ 0.50 in the low spin state and is larger than 50% from (3 ~ 0.20 to ~ 0.70. As is shown 
in the previous subsection, the intrinsic wave function on the J 7r =0 + curve starts to change 
its structure from a deformed mean- field-like structure to a cluster structure around /5=0.5. 
Therefore, unlike the K w =0i band, the K T =0f band is a mixture of the cluster structure 
and the deformed mean-field structure. Indeed, the wave functions of the ground band have 
smaller overlap with the a+ 16 model space than the K n =0i band fTablelLT]). This tran- 
sient nature of the ground band is reflected in the expectation value of the spin-orbit force, 
which amounts to about -5.0 MeV in the ground state and increases as the angular mo- 
mentum increases. The increase of the spin-orbit force contribution is a consequence of the 
structure change along the yrast line or the anti-stretching phenomenon, which is discussed 
later. 

The transient nature of the ground band affects the intra-band E2 transition probabilities 
TablelTV|). They are successfully described quantitatively by the cluster models (RGM 
and coupled channel OCM [3|), but they need to introduce a small effective charge. 
When the model space is expanded to include the a+ 16 O(0 + , 3~, 1~) channels (multicluster 
GCM) 22(, the deformation of the mean-field is described to some extent and the result is 
improved. In the present calculation, the E2 transition probabilities of the ground band are 
described well without any effective charges. Therefore, we can conclude that the general 
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trend of the intra-band E2 transition probabilities are well described qualitatively by the 
a+ 16 cluster component in the wave functions and the existence of the deformed mean- 
field component (or the cluster dissociated component) slightly enlarges the intra-band E2 
transition probabilities. 

Another important feature of the a+ 16 cluster structure is the anti-stretching phe- 
nomenon, which is studied by many cluster models Q. As the angular momentum in- 
creases, the a+ 16 cluster structure moderates and the average distance between a and 16 
becomes smaller. Indeed, in the present result, the position of the maximum peak of the 
ct-RWA moves slightly toward the inside of the nucleus as the angular momentum increases, 
though the tendency is rather mild. However, we wish to emphasize the difference of the 
anti-stretching phenomenon described by the a+ 16 cluster models and the present model. 
In the a+ 16 cluster models, the anti-stretching phenomenon is described by the decreasing 
distance between a and 16 O. However, in the present model, as the angular momentum 
increases, the a+ 16 component diminishes and the spin-orbit force acts strongly. In par- 
ticular, the 8^ state is oblately deformed and four valence nucleons are aligned, which leads 
to a large contribution of the spin-orbit force which amounts to about -10 MeV, though 
a small amount of the a+ 16 component survives outside of the nucleus. We also note 
that this remaining small amount of a+ 16 component outside of the nucleus leads to the 
relatively large a decay widths of the J 7r =6^ and 8+ states (Table. UJ). 

3. K w = 0+ band 

The K 7r =Q~l band has a prominent a+ 16 cluster structure in which the relative motion 
between clusters is 2ftw excited (N = 2n + L = 10), while a and 16 clusters are not 
excited. In the present calculation, the obtained second positive-parity band members have 
large amounts of a+ 16 component and large a decay widths. Also, the relative motion 
between a and 16 has the principal quantum number N = 10. Therefore, we identify 
this band as the K*=Qi± band. The K n =02 and Og" bands which do not have the a+ 16 
structure are not obtained in the present calculation, but we consider that we will obtain 
these two bands when we make the variational calculation and the GCM calculation in 
the model space which is orthogonal to the K n =0 + (J 7r =0 + , 2 + , ...) curves obtained in the 
present calculation. The calculated and observed excitation energy of the O4" state are 8.44 
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FIG. 8: The GCM amplitudes for Of (solid line) and 0^J~ (dashed line and dotted line) states. The 
definition of the GCM amplitude is given in the text. 

MeV and 8.3 MeV, respectively. 

We consider it worthwhile to discuss how the prominent a+ 16 cluster structure of this 
band compared to the ground band is derived in the present calculation. As explained in 
section |HJ we include all of the obtained intrinsic wave functions as the GCM basis. The 
prominent a+ 16 cluster structure is due to the considerable contribution of the 3>j nt (_)(/3) 
which was obtained by the variational calculation after the projection to negative-parity. 
When we include only P^ K=0 ^ int ( + )((3) in the GCM basis, the content of a+ 16 components 
of K n =0^ band members is around 40~70% (bracketed numbers in the third column of 
Table |HJ and comparable with that of the ground band. In contrast, when we include 
the P^ K=0 ^int(-)((3) as well as P/ f H ^ =0 $j rit ( + )( / 3), the a+ 16 cluster structure of this band 
becomes prominent (numbers in the third column of Table lH|) . In the outside of nucleus, the 
contribution from P^ K=0 ^ int ^)((3) (dotted line in Fig. |HJ) becomes larger, while in the inside 
PMK=Q$int(+)(P) (dashed line in Fig. (HJ) is dominant. Therefore, in the outside of nucleus, 
the a+ 16 cluster structure is prominent, which leads to a large a decay width. However, 
in the inside, the mean-field-like structure contributes to some extent and it leads to the 
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TABLE V: Observed and calculated intra-band E2 transition probabilities B(E2; Jf — ► JJ ) |e 2 fm 4 1 
within the K n =2~ band. For comparison, the results of the j-j shell model (SM) |2l| and 



(a+ 16 0)+( 8 Be+ 12 C) coupled channel OCM (OCM) 



2Cj are also shown. The effective charges 



which were used in each calculation are given in the bottom line. 



T w — > 


TTT 
J f 




Olvl 




AMD-l-CtCIM 

ii 1VX 1J ^ vi V_y ivx 


3f - 


2~ 


113 ±29 


97 


108 


102.8 


4~ 


sr 


77 ± 16 


75 


77 


77.8 


4" 


2~ 


34 ±6 


36 


34 


38.5 




4" 


< 808 


44 


45 


84.5 




sr 


84 ± 19 


48 


49 


56.6 


6~ 


sr 


32 ± 13 


32 


34 


29.9 


6- -> 


4- 


oo_ 13 


51 


67 


64.0 


Se/ 


e 




0.8 


0.069 






existence of the small non-a+ 16 component. It is interesting that though PMK=o^int{-){f3) 
is less bound than P£r K=0 § int ( + )(j3), PMK=o®int(-)(P) contributes to this band. 



4. K n = 2~ band 

The K n =2~ band is the lowest negative-parity band. It starts from the 2~ state at 4.97 
MeV and is observed up to the 9~ state, and exhibits a good rotational spectrum. The 
obtained excitation energy of the 2~ state is 6.51 MeV and we have calculated up to the 6~ 
state. 

This band has a quite different character from the other three bands. From the shell 
model calculations, it is known that this band has a (lp)~ 1 (sd) 5 structure with respect to 
the 16 core and the 57/3(8,2) configuration is dominant, which means the absence of the 
a+ 16 cluster structure and a large contribution from the spin-orbit force. Indeed, in the 
present result, the GCM wave functions of this band members consist of the wave functions 
on the K n =2~ curve and in the natural parity states (n = (— ) J ), the mixing with the wave 
functions on the ^"=0" curve is quite small. As discussed in the previous subsection, the 
wave functions on the K n =2~ curve have a deformed mean-field- like character. Therefore, 
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the spherical-basis AMD gives much higher excitation energy of the 2~ state, 11.3 MeV. 
Though the character of the band is quite different from that of the ^^=0^ band, the 
B(E2) transition probabilities are again described well without any effective charges. 

V. SUMMARY 

In this study, we have presented the details of the new theoretical framework of deformed- 
basis AMD. This framework enables us to describe the coexistence, competition and mixing 
of the deformed mean-field structure and the cluster structure and gives us the proper 
evaluation of the spin-orbit force. 

We have applied this framework to the 20 Ne nucleus. On the analysis of the obtained 
energy curves, three different characteristics of the nuclear structure appeared. The J n =l~ 
(K n =0~) curve has an almost pure a+ 16 cluster structure, in which the spherical single- 
particle wave packets are localized into two parts. In contrast, the J 7r =2~ (A 7r =2~) curve 
has a deformed mean-field- like character in which the significantly deformed single particle 
wave packets are gathered around the center of the nucleus. The J 7T =0 + (A' 7r =0 + ) curve has 
a mixed character. As the nuclear deformation becomes larger, the wave function gradually 
changes its structure from a deformed mean-field-like structure to a cluster structure. 

By superposing the obtained wave functions, four rotational bands are obtained in the 
low-energy region. The K n =0i band has an almost pure a+ 16 cluster structure. Observed 
large a decay widths of this band members are reproduced and the obtained a-RWAs show 
a quite similar trend to those obtained by other cluster models. These facts mean that the 
cluster structure can exist as an independent degree of freedom of the nuclear excitation, even 
though the deformation of the mean-field and the effects of the spin-orbit force are included 
in the model space. The ground band has a mixed character of the deformed mean-field 
structure and the cluster structure. Because of the mixture of the deformed mean-field 
character, the obtained intra-band E2 transition probabilities are enlarged compared to the 
cluster models and the observed values are reproduced without any effective charge. The 
structure change along the yrast line is also investigated. As the angular momentum becomes 
larger, the deformed mean-field component in the wave function increases and the spin-orbit 
force acts strongly. By including the wave functions which have an almost pure a+ 16 
structure as the GCM basis, the K^=0^ band acquires a prominent cluster structure and 
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observed large a decay widths are reproduced. The K w =2~ band has a different character 
from the other three bands. The wave functions of this band member have a deformed mean- 
field character. The intra-band E2 transition probabilities are again reproduced without any 
effective charge. 

Though these four bands have different characteristics, deformed-basis AMD successfully 
describes the coexistence and mixture of the cluster structure and mean-field structure in the 
same nucleus. Therefore, we believe that deformed-basis AMD is one of the most powerful 
approaches to describe the interplay of the mean-field structure and the cluster structure 
and to study the cluster structure and exotic cluster-like structure in medium-heavy and 
neutron-rich nuclei. 
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APPENDIX: REDUCED WIDTH AMPLITUDE OF THE AMD WAVE FUNC- 
TION 



The calculational methods for the a reduced width amplitude (a-RWA), a decay width 
and the square amplitude of the a+ 16 component are briefly given. The method for 
evaluating these quantities from the AMD wave function was developed by Y. Kanada- 



En'yo 



and the reader is directed to references 



for more details. 
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The a+ 16 system is generally expressed by the RGM-type wave function. 

S£«o = ^{xj(r)lj O (f)0(a)0( 16 O)}, (A.l) 
n = v/16!-4!/20!, tt = (-) J . (A.2) 

Here, A is the antisymmetrizer, r is the relative coordinate between a and 16 0, <fr(a) and 
0( 16 O) are the normalized internal wave functions of the clusters. Xj{ r ) is t ne radial wave 
function of the relative motion between a and 16 0, and so normalized that <3> J . 16r . is nor- 
malized to unity. 

The normalized deformed-base AMD+GCM wave function of 20 Ne is divided into 
the a+ 16 component and the residual part . 



*m = a^ +16o + VT^^. (A.3) 

is also normalized and orthogonal to the $^ie . We introduce the projection operator 
Vl which projects out the a+ 16 component from the 

= a n A{xj(r)Yjo(f)<p(a)<P( w O)}. (A.4) 

The practical formula for Vl used in this study is given later. Using this projection operator, 
the squared amplitude of the a+ 16 component Wj of $^ is written as, 

Wj = \a\ 2 = (*ft\P L \* J J)- ( A - 5 ) 
The a-RWA yj(a) is defined as 

yj(a) = -(^W(a)^ 6 0)|^) 
no r z 

= -(^^F JO (f)0(a)0( 16 O)| 

x|*£«o>> (A.6) 

and its squared amplitude |a3^(a)| 2 is the probability to find a and 16 clusters at the inter- 
cluster distance r = a. If we expand Xj{ r ) with the radial wave function of the H.O. with 
the width parameter 7, R^j(r,j) as 

Xj(r) = e NjRNj(r, 7), (A.7) 
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the a-RWA is written as 

y.j{a) = a n N e NJ R N j(a, 7). (A. 8) 

N 

This means that when is obtained, y(a) is expanded by the H.O. wave function. Here, 
the eigenvalue of the RGM norm kernel /i^ is defined as 

Hn= (Rnj(r)Y J0 (r)cf>(a)<f>( le O)\ 

^{i2nj(r)y JO (f)0(a)0( 16 O)}), (A.9) 

and its value is analytically evaluated. We calculate a decay width T a using the a-RWA. 

T Q = 2P/(a) 7 2 (a), 
PAa) 



Fl(ka) + GKkaY 

7 2 (a) = ^l«^(a)| 2 , (A.10) 

where Fl and Gl are the regular and irregular Coulomb functions, k is the wave number of 
the resonance energy, and a is the channel radius, which is chosen to be a = 6fm. 

To evaluate Wj and Xr( a )> we have approximated the projection operator Vl with the 
set of the orthonormalized angular-momentum-projected Brink wave functions. We start 
from the Brink wave function (fisiRi) i n which a and 16 are represented by the S^-limit 
wave functions which are located at the points (0,0, ^§-R«) and (0,0, — ^Ri), respectively, 

16 4 
<p B (Ri) = n A{tp(a,—IU)M*O,-—1L i )}, (A.ll) 

Ri = (0,0,^). (A. 12) 

The width parameters of a and 16 O are taken to be the same value v for simplicity. By 
separating the center-of-mass wave function w(Xg), the internal wave function (pc{Ri) is 



written as 



<p B {Ri) = w(X G )ip c {Ri), (A.13) 

<pc(Ri) = n o ^{r(r,R i , 7 )0(«)0( 16 O)}, (A.14) 

r(r,R,, 7 ) = f^j e"^) 2 , (A.15) 

lu(X g ) = f-^J exp(-20z/X^), (A.16) 
1 20 Ifi 4 

i=l 
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where T(r, Ri, 7) is the relative wave function between a and 16 O of the Brink wave function. 
The angular-momentum-projected Brink wave function (f) J c (Ri) is obtained from the <j)c{Ri), 

ip J c {Ri)=qiP^c{Ri). (A.18) 

Pq Q is the angular momentum projector and qf is the normalization factor. The orthonor- 
malized set of the angular-momentum-projected Brink wave functions (p J a is described by 
the linear combination of the <Pc(Ri), 



and the coefficients p a and Q a are the eigenvalues and eigenvectors of the overlap matrix 

B ij = (<p&R i )\<pl } {R j )), 

BijCi a = p a Cia- (A. 20) 

If enough numbers of the basis wave function (p J a are used, the projection operator Vl is 
approximated by the set of tp J a) 

a. 

= J2 B 7'\v J c(R*))(v J c(Ri)\- (A.21) 

ij 

In the present calculation, we have employed 22 Brink-wave functions in which Ri is taken 
as Ri = 0.5, 1.0, 1.5, ...,11.0 fm. 

Using Eq. (TOTJ) . Eq. (|A~4|) and Eq. (jA~5j) are rewritten as 

Vl$% = an o ^{xj(r)lj O (f)0(a)0( 16 O)}, 

ij 
ij 

x(^\cp J c (Ri)). (A.23) 

The expansion coefficient of Xj{ r ) an d y( a ) is obtained from Eq. (jA.22|) . We expand 
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the wave function of the relative motion 7(1-, Rj, 7) of Eq. (| A. 24(1 



r(r,R i , 7 )=e-^ 2 ^ 



NJm 



2^j;Ajy Jm (4)^j(r,T)ijm(r)j, (A.24) 

Aj ~ { j V (N — J)!! • (N + J + 1)!!' (A ' 25) 
And using Eq. (IA.14jl and ()A.24jl . Eq. (IA.18jl is rewritten as 

ip&Ri) = n A{ X f(r)Y JQ (r)<p(aW 6 O)}, (A.26) 

Xf(r) = qfe-^ £ W^ A »R NJ {r). (A.27) 

AT V -/V. 

By inserting Eq. (|A.26|) into Eq. (|A.22jl . we obtain x.j{ r ) as a superposition of Xj ( r )> 

^ -E^MIO^M^W- (A.28) 



Inserting Eq. (IA.8j) and (IA.27j) into Eq. (IA.28jl . e^v is given as 



xg/e -i«? (7^1 . (A 29) 
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